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Under certain conditions the number of photons radiated classically by a charged particle follow- 
ing a prescribed trajectory can be finite. An interesting formula for this number is presented and 
discussed. 

Ctf ■ 

I. FINITE NUMBER OF PHOTONS AND THEIR FORMULA 

o ■ 

• In the classical theory of radiation by a charged particle one calculates the energy radiated into the electromagnetic 
field. Indeed from the purely classical point of view the energy is practically the only quantity there is to calculate. 

i-H However, in certain problems one may come upon the idea of finding the number of photons radiated by a charged 
^ particle following a given trajectory. Although the photon is a quantum concept and so the question of finding the 
, number might be thought to involve quantum mechanics, it does so only in the most minimal way. Given the energy 
radiated into a given mode of the field, it is only necessary to use Planck's relation and divide by the frequency u> to 
find the number n; thus the problem remains an essentially classical one. 

It might be objected that n can be infinite while the energy is finite, as in the well-known "infrared catastrophe" . 
£N| ' True, but as we shall see, there is an interesting class of cases where this is not the case. In particular when the 

• trajectory of a charged particle begins and ends with the same vector velocity v, n is generally finite. This is because 
' the "infrared catastrophe" results from the difference in the long flight paths for the in-and-out going particles, and 

> , when they are the same, the "catastrophe" is averted. 

On the other hand, bremsstrahlung calculations often have an "ultraviolet" or high frequency divergence, due to 
D | the sudden appearance or deflection of a charge [Q . Evidently one cannot emit an infinite number of finite energy 
' photons, so this divergence must be an artifact of the calculation. This difficulty can have either a quantum or classical 
resolution. In the quantum solution, as in the feynman graph technique, one takes into account the energy-momentum 
conservation usually neglected at the classical level, and a suppression results for high energy photons. However, and 
more of interest to us here, there can also be a classical resolution: the ultraviolet divergence results from abrupt 
changes in the trajectory, and if the velocity changes or accelerations are sufficiently smooth, there is no divergence. 

Therefore we have the interesting situation that for a smooth trajectory, beginning and ending with the same veloc- 
ities, the number of photons n radiated according to a classical calculation should be finite. But n is a dimensionless, 
lorentz-invariant quantity. Thus there ought to be some simple formula for it, a relation mapping the path of a particle 
in space-time to the real number, n. 
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This relation is: 



a f f , 1 



n = — 



dx^-^dx' (1) 



IT J J —»S L 

In this formula x and x' are four-dimensional coordinates, referring to points on the space-time path of the charged 
particle, so that dx^ is a 4-vectorial element of the path (Fig 1). One may introduce the proper time or invariant 
path length r and the four- velocity it M = dx^/dr to also write the expression as 



- = j J i^i? iTiT ' =-«!! "» <T ' ) **' <2) 

Si e is the four-distance between the points x, x' in the following way 

^ = (<-<' + z £ ) 2 -(x-x') 2 (3) 
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As we shall explain shortly the it is necessary to make sense of and to properly define Eq jjj . 

It is interesting that the Sommerfeld fine structure constant a = f- ~ 1/137, involving Planck's quantum constant, 
appears. It arises, however, not as a coupling constant, but rather from the conversion of the classical energy to 
quanta, where we finally must divide by Tilj. The expression is a non-perturbative. If the particle is multiply charged 
with Q electron charges, the formula should be multiplied by Q 2 . Observe that by interchanging t and t' in the 
integrations one can reverse the sign of it to show that n is real. The expression is not an integer since it represents 
the average or expectation value of the number of particles radiated. 

t 
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FIG. 1. Path in space-time 



II. DERIVATION 



We can derive Eq |l| from the standard treatments of classical radiation theory ||, j|] where one calculates j(k), 
the fourier components of the three-vector current, and then finds the energy radiation rate ~ |j| 2 sm 2 #, where 9 is 
the angle between k and j. Dividing by u (units h = 1) gives n. 

To obtain the total number, we sum over all modes k of the radiation field, 

sm 2 0|j(k)| 2 - d 3 k . (4) 

One could try to perform the integrations for every particular trajectory, but we would like to find a general formula 
in terms of the path itself. Thus we attempt to perform the d 3 k integration first, before the fourier transform. This 
gives at first a not well-defined expression, which leads to the need to introduce the it. 

To see how this comes about, we observe that for smooth paths x(f) the fourier transform j (k) is strongly convergent 
at large u>. When we introduce the current density for our classical point particle as 

j(t,x) =ev(i)£ 3 (x-x(t)) (5) 

we have j(k) = e J dtve 1 ^*^^^ for the fourier transform. Since the path of the particle is time-like the exponent 
can never become zero, and if x or v = ^-x have no jumps or kinks the integral vanishes rapidly for large lu. Eq Q] 
can be rewritten in a more covariant way using sin 2 9 = 1 — cos 2 9 and noting the current conservation relation 
ujj = ujcos9\j\, so that Eq [Q is essentially the integral of the four- vector squared j 2 — j 2 = —j 2 . Using dt v = dx 
one has a compact line integral expression for , namely 

j M (k) - e [ dx^ kx . (6) 
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where kx is the four dimensional scalar product k v x v = hot — kx(i), with fco = u> = |k| 

Squaring Eq || and writing in the constants || we get for the integral of —j^(k) over all modes of the radiation 
field: 

n =~£^J dZk l J J dx^dx'^-^ , (7) 

Firm in the belief that we actually have a well defined, convergent expression, we allow ourselves to do the d 3 k integral 
first and to use some ie manipulations in dealing with seemingly ill-defined expressions. 

After a few steps we get to the expression ujduje luJ ^ t ~ t )~ Rcos0 ) ; where R = |x — x'| and where here 8 is the 
angle between k and the vector (x — x'). Now, since we believe that our expressions are rapidly vanishing at large 
uj and lu is always positive, it shouldn't hurt to add an ie in the exponent Q, to make the replacement (t — t') — > 
(t — t' + ie). The integral is then convergent, and integrating over dcosQ and keeping track of the ie, this leads to 
J^j.lgifx-x') = _4tt^- and so Eq §. 

Although these manipulations are similar to those with the invariant D functions in quantum field theory, and our 
function 1/S 2 resembles a D function, our assumption of smooth paths has lead us to an "ie prescription" in position 
rather than in momentum space. 



III. ELIMINATION OF e 



The difficult job now becomes the evaluation of the ie in Eq Q. We look for guidance to the theory of generalized 
functions || where one has the relation 

+oc m dt f 00 ittiA dt (8) 
co ww ~ J-oo ~^^ dt ' (8) 

where we have specialized to the case of an even function f(t) = f(—t). 
One way of understanding this relation is to note that, for e non-zero, 

which just follows from explicit integration. Thus we have simply subtracted zero in Eq and have chosen the 
coefficient of this zero term in such a way as to cancel the singularity of the integrand. With f even, f(t) — f(0) ~ t 2 , 
and the expression is indeed finite at the singularity. 

Can we apply the same idea here? That is, by examining the neighborhood of the singularity for small but non-zero 
e we see that our expression is finite and e independent. We would thus like to subtract something which is zero for 
finite e and which will regulate the singularity at t = t' in Eq Q, enabling us to finally dispense with e altogether. The 
problem, however, would seem to be much more difficult than in Eq H. There we simply had to adjust one constant, 
namely /(0). Here we need a function f(t,t'), such that when it is integrated over l/S 2 (t,t') gives zero and then can 
be adjusted to cancel the numerator dx^dx'^ when we make the replacement dx^dx'^ — > dx^dx' (1 — /). Furthermore 
it should do this for any path x(t) we care to put in the formula. 

This last requirement rules out, for example, the at first likely- looking candidate replacement dx^dx'^ = 
u^u'^drdT 1 — > {u^u 1 ^ — l)drdr'. Since u^u'^ — > 1 for r — > r 1 this would regulate the singularity and looks promising. 
But it doesn't work, J J drdr' 1/S 2 (t, t') is obviously path dependent and cannot be zero for all paths. 

This, however, suggests a solution. If we can find something for / which is a total derivative and so can be integrated, 
it would depend only on the end points of the path and not on the path itself. To this end, note the following useful 
relation. Let G be some function of S 2 (t,t') = where A p = x^(t) ~ 2^ t (Y'), then 

d T d T> G{S 2 ) = -G"4(A^)(A^) - 2G'^^ , (10) 

where G' and G" refer to first and second derivatives with respect to S 2 . 

Since the expression is a total derivative and its integral is independent of the path and is to be interpreted as zero, 
we may look for some choice for G that leads to ~ 1/S 2 for r — > r' and so might cancel the singularity. This occurs 
in fact if we consider G = InS 2 . The above equation becomes 



3 



«A, MS V + >,£)(A.£)-£££. (ID 

Since ^r- — > for r — > r' the expression ~ 1/5 2 at the singularity. Now we want to add this to Eq that is to 

/ -^r^p-(l/ S 2 )dTdr' , with some weight such that the singularity at r = t' is cancelled. To determine this weight 

note that since the four velocity satisfies (^r) 2 = 1 the rhs of Eq goes to +2/S 2 for r — > r'. Hence we must 
subtract 1/2 of Eq |T| from Eq g to remove the singularity. Doing this we end up with the following nice relation. 



drdr' = 2 I I 7 M "^ 0^ -3^ drdr ' ( 12 ) 



dr S£ rfr' dTdT ~\ I dr S 2 dr 



or alternatively 



1 , / n f f , *W s' 2 



dx^K = 2 J J dx » S 2 d < ( 13 ) 

The notation is meant to indicate that while the ie is in S 2 on the left, it is not needed on the right. 

We can rewrite these equations in an interesting way by noting that is like an average 4-velocity connecting 

the points r, r', which we might call in analogy to the instantaneous 4-velocity = 

^(r,r , )4 ^" (T) s I " (T ' ) (14) 



Then we can write 



o 



U ( T ) ■ <^')-{U ■ u(t)) (U ■ u{t')) 



■■ r 2 J J u » s 2 Uv dTdT ' = k 2 J J — s 2 drdr ' (15) 

We can now introduce the notion of a "transverse vector" u T associated with the points x,x' at r, r'. This is the 
vector u with the "longitudinal" part, that is the component along U, removed: 

uI(t,t')=u im (t)-U im (U-u(t)) uI(t\t)=u i1 (t')-U i1 {U-u(t')) (16) 

With these definitions, u^(x, x')u^(x', x) = u il {x)(5^ lJ — U l JJ v )u v {x'), which is the numerator in Eq |l5| and we can 
write the above equations in terms of products of "transverse vectors" : 

a 2 f ful(r 7 r')^ul(r',r)drdr'. (17) 
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Hence we can say that our integral for n represents the "interaction" of pairs of transverse vectors along the path of 
the charge. The possible singularity at S — is absent because u t (t,t') vanishes for two points very close together. 



IV. PROPERTIES OF THE QUANTITIES 

We summarize some properties of these quantities: U and u are unit vectors, u 2 = U 2 = 1. From their definitions 
U and u T are orthogonal v/^(t,t')U^ = u^(t',t)U^ = 0. Note for straight line motion or more generally for any 
smooth path as r — > r' 

U» -» u„ (18) 

Unlike u and U, u T is a space-like (or zero) vector. This follows from the fact that there is a frame where the time 
component of u T vanishes, namely the restframe of the time-like U . 

We note a point concerning the definition of U in Eq [Q. This point is irrelevant in all those expressions where 
S or U appears quadratically, but we mention it for consistency. Namely, we would like U to resemble the velocity. 
Therefore it must be understood that S(t,t') = \f~S 2 is "directed"; an odd function, positive for r > r' and negative 
for r < t', like (r — r'). It implies a definition of U such that U(t,t') — U(t',t) and that U is "forward pointing", 
i.e. Uq is always positive. 
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An interesting expression for u T follows from the definition of U in Eq |Q, reflecting the fact that the derivative 
of a unit vector is transverse to itself: 

i uI(t,t>) = drU^r') - \ <(r',r) = d T ^(r,r') (19) 



With this, we can also write Eq |17| in some different looking ways: 



// - -- 2 / / d T U^{T,T') dr.U^T') drdr' (20) 



Or, since U 2 = 1 and so d T d T ' (V ]JJ M ) = we can also write 



2 / / r.Ar.r) ,L ( L, f - /( , T . T > ) drdr' (21) 

Another variant results if we note that Eq |2(J looks like the cross terms of a quadratic expression: 2d T Ud T >U = 
\[{d T U + d T 'U) 2 — (d T U — d T 'U) 2 ]. Introducing the sum and difference variables r + = i(r + r') and r_ = i(r' — r) 

n=^ I I [u> T r.r - «) 7 ? ■ r] ,/r + </~ . !22i 

where because of the rotation of the coordinates in the r, r' plane, the integration boundaries now have the form of 
a diamond instead of a square. 



V. INFRARED BEHAVIOR 

The ie takes care of the high frequency or ultraviolet behavior, but there might be questions concerning the infrared 
region. We see this very simply if we consider the stationary particle with v = 0, for which n of course should be 
zero. Eqjl} leads to 

dtdt'- - — (23) 

There are two kinds of limits implied here e — > 0, and some upper/lower limit of the integration L — > oo, and the 
integral can depend on how the limits are taken. If we simply apply Eq || we get of course zero, as desired. On the 
other hand if we interpret the limits at ±oo by, say, integrating over some test function which is constant up to some 
large number L and then drops off, we can get an answer involving eL, which depends on the order in which we take 
e — * 0, L — > oo. The problem originates in the fact that we are turning off the charge at some L, violating charge 
conservation and in the process producing some photons. We are thus to handle this limit by remembering that e is 
to be kept finite to the very end, and so we first mean eL — > oo. Alternatively one can subtract the integral for the 
stationary particle, as in Eq |2^] to regularize the infrared behavior. 

Naturally if we use the form without the ie, say from Eq |l5| we immediately get zero for the stationary particle. 
More generally, we can investigate the contribution of the final and initial long straight paths for a moving particle in 
Eq The dan gerous regions for the infrared behavior are large r, r' as the particle comes from or goes to infinity. 
For both r's very large or both very small, we have in view of Eq [Q and u 2 = 1 that the numerator becomes 
u 2 — u 2 u 2 = 0. Similarly if one r is at very early times and the other a very late times and the velocities at these 
times are the same, there is the analogy to Eq |18f 

UJr,r')= Ufl + 0(^— ) (24) 



where T is the finite time period where the charge was not in uniform motion. For large (r — t'), U — > u and the 
numerator will again tend to zero, as long as the initial and final velocities are the same. It is interesting that the 
elimination of both the long and short distance singularities can in a sense be attributed to the same relation, Eq |L8| . 
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VI. NON-RELATIVISTIC LIMIT 



Consider that class of paths whose tangents are roughly parallel to some straight line on Fig (1), meaning that the 
ordinary velocity v is always close to some typical or average velocity. By making a lorentz transformation (under 
which n is invariant) so that this typical velocity is zero, we have as a first approximation in the ordinary velocity the 
non-relativistic limit, leading to expressions quadratic in velocities. 

We could make a stab at the non-relativistic limit by taking our basic expression Eq and just naively expanding 
it for small velocities. We introduce V, the average velocity vector connecting two points on the curve 

V= *(*)-*(Q. (25) 

t-t' y J 

V is a symmetric quantity V(t,t') — V(t',i), in parallel with our earlier definition of U and becomes equal to v 
when t — > t' . V plays a role analogous to U except that it does not have a fixed length and so V 2 has non-zero time 
derivatives. Expanding in Eq [jjj we can try to write 

1/Sl « {t _ t ! +ie) 2 (! + + V *W) + -) ( 26 ) 

where V ie = ^*^*2=) so that 



n = 



I J J dx^dx'^^ J J M'(l-v(t)v(0) / + (1 + V*(t,t')+V*(t,t') + ...) ■ (27) 



Now, using Eq M the "1" term vanishes and we have 



ol f f , , , 1 



n = — 



I J dtdt'- — ^ [V 2 - v(t)v(t')] (non- rel. limit) (28) 

We have dropped the ie since the expression is non-singular with V — > v as t — > t 1 . 

This is a not unreasonable-looking expression. Indeed, if we expand it for t w t' we get 

[F 2 -v(t)v(t')]« \a 2 (t-t') 2 -, (29) 

showing that the singularity is cancelled and that the leading terms of the expression are positive and exhibit the 
familiar connection between acceleration squared and radiation. The (...) includes terms involving time derivatives 
of the acceleration a (a = 4|) as well as higher order terms in (t — t') 2 . 

One might have some qualms about the carefree ie manipulations, and we indicate how to arrive at Eq p8| by 
straightforward application of our more conventional formulas. 

First note that U can be written in terms of V in the usual way relating a three-velocity and a four-velocity, 
U Q = l/Vl - V 2 , U = V/a/I - V 2 with the non-relativistic limits U « 1 - 1/2V 2 , U w V. we find from either 
Eq Q or Eq Q to leading order in V 

n=-2j J dtdt' d t V(t,t') d t >V(t,t') (non - rel. limit) (30) 

Now in analogy to Eq we have from the definition of V 

ftvco-^ a ,v<u<) = -^ (31) 

and Eq [|o| becomes 

n=-2 I |«'[-v(()v((') + (v(()+v(('))V-y 2 ]- — ^—^ (non -rel. limit) (32) 

7T j j (t-t'r 

One could be satisfied with this formula as it is, but to bring it into the perhaps simpler form Eq pjj] , note the 
following identity: \d t d t >V 2 = [-v(t)v(t') + 2(v(i) + v(t'))V - 3F 2 ](< - f') -2 ) wnich follows from differentiating 
V 2 = (x-x') 2 (t-i')~ 2 twice. We split this into two parts: \d t d t >V 2 = [-v(t)v(t') + (v(t) + v(t'))V- F 2 ] (t- t'Y 2 + 
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[(v(t) + v(i'))V) - 2V 2 ](t - t'y 2 , where we make the split so that the first part corresponds to Eq Q. Now d t d t 'V 2 
is a total derivative whose integral may be set to zero. Therefore the integral of the first and second parts represents 
the same quantity with opposite signs and we can write: 



dtdt'[-w{t)w{t') + (v(t) + v(t'))V - V 2 



1 



(t-tf) 



A2 



dfdt'[(v(<) +v(t'))V) - 2V 2 



1 



(*-*') 



M2 



(33) 



Taking the one-half the sum of the two forms we finally obtain, after much labor, Eq pSy . The cavalier ie manipu- 
lations were certainly a lot quicker! 

Despite the familiar acceleration squared in Eq [^9|, we shouldn't expect that n can be represented simply by an 
integral of some local quantity along the path. There is the (...), and all our expressions are bi-local in the time. 
The photon is a non-local concept and a certain space-time interval is necessary to define it. This corresponds to the 
distinctive property of relativistic local field theory that while one has local constructions for quantities like charge 
density or energy density, there is in fact no local quantity for particle number or photon density. Indeed, the need 
to have some space-time interval to define a particle, leads to the concept of the "formation zone" || , which can be 
used to understand certain phenomena like the absence of "cascading" for particle production on nuclear targets. 



VII. SIMPLE CASES 



With the simple Eq |2^] in hand we can proceed to calculate a couple of concrete examples. 

Dipole radiation: We first take the classic problem of dipole radiation. Let a charged particle be oscillatin g in 
one dimension according to x = Xosinflt, so v = x^ftcosftt. Changing variables fit — > t and similarly for t', Eq |2q| 
becomes 

O- , f f , , / risint — sint') 2 cost cost' , . . 

n=-(x flfj Jdtdt' ^^^.^ - {t _ t/ + ie)2 ] (34) 

Although the ie is not necessary since the combination of the two terms gives something non-singular, it is convenient 
to keep it since it allows us to handle each term separately. Carrying out, say, the t integral first and using relations 
of the type / dt u™, 8 ^ = 27re~ lt ' or / dt {t J™ + \ e)i = TT{i/Q)e~ u ' and / dt jr^r^p = 0, we arrive at 

n=^-(x fl) 2 J dt'[^-sint' - (-Tr)cost']e- lt ' = a(x Q fl) 2 ™ (35) 

where T is the length of time the particle is in motion. The number of photons generated increases linearly with time, 
as was to be expected. We have neglected a contribution, not proportional to T, connected with turning the motion on 
and off. If we now multiply by ft to find the energy and divide by T to get the power, we obtain Power = a^(xofl) 2 fl 2 , 
which is the classical formula for dipole radiation averaged over a cycle j| , || . Not surprisingly we recover the classical 
result, as was our starting point. 

Note however, that in general the energy radiated and the number of photons do not stand in direct relation since 
the oscillating charge produces higher harmonics in addition to the fundamental at frequency ft jjj. However, in the 
nonrelativistic limit v/c — ► these higher harmonics become negligible (since they are a retardation effect), and we 
expect the energy radiation to be simply proportional to n. 

This example shows that despite our requirement that the initial and final velocities be equal, the method need not 
be of purely academic interest for practical calculations. If the effects connected with turning the motion on and off 
are negligible compared to some main effect, we can always return the particle to, say, zero velocity, while retaining 
the main effect. 

Smooth deflection: Instead of an oscillator which is on for a long time we can consider a charge undergoing a 
smooth deflection, for example x — in . , ,},. ^ so that v — —2^-^- . , .},. vi . This leads to 

' v u l+(t/to) J to to l+(t/to) 
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VIII. FURTHER QUESTIONS 



Eq |l| should have some general symmetry properties with respect to changing the path. There are the evident 
invariances under lorentz transformations, translation, 3-D rotation, reflection, and time-reversal. Since there are no 
dimensional quantities except the path itself involved, there is also an invariance under rescaling of all 4-coordinates 
simultaneously, x M — > Xx^. That is to say, if the path is expanded in space and time proportionally so that the 
velocities remain unchanged, n is unchanged, as we see in the examples. It would be interesting to know if there are 
further invariances and what the full invariance group is. 

Also, we might consider the problem for gravitons instead of photons. Presumably one will find, in analogy to Eq 
and in view of the tensorial character of the source of gravitons, n ~ Gm 2 J J u A1 (T)w^(T)-^-u AI (r')/iy(r') drdr' where 

G is the gravitational constant and in the mass of the radiating particle (recall h, c = 1, so Gm 2 — (m/M p i) 2 ), but it 
would be interesting to investigate this more closely. 

Finally, a number of interesting mathematical problems suggest themselves. Our n gives an invariant characteriza- 
tion of the "wiggly-ness" of a curve. There appears to be no reason why it should not be also used in euclidean space 
where S 2 (t, t') = T,{ Xi {t) - Xi (t')) . The main difference would seem to be that the curve can now go "backwards" , 
opening the possibility of closed curves. As well, there is the possibility of new singularities when two parts of the 
curve, with remote values of r, come close together. 

For a closed plane curve it is plausible that the minimum value of n obtains for the circle. Then defining n — 
—2 J J ^2U T u' T drdr' ( we use Eq |TJ, leave away ^, and the natural sign is now minus), we thus conjecture that 
the minimum value of n for any plane curve is 2tt 2 , which is what we obtain for the circle. 

Also, there may be a "topological" aspect to n, connected with knots. For an open path the minimum value of n, 
namely zero, is reached for a straight line. If the path has a knot, however, it cannot be continuously deformed to a 
straight line, and must go "backwards" somewhere, suggesting that the minimum value of n continuously attainable 
is related to the presence of knots. 

I am grateful to Y. Frishman, D.Maison, and V. Zakharov for many and helpful discussions. 
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